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Abstract 

The structure of hierarchical networks in biological and 
physical systems has long been characterized using the 
Horton-Strahler ordering scheme. The scheme assigns an in- 
teger order to each edge in the network based on the topology 
of branching such that the order increases from distal parts 
of the network (e.g., mountain streams or capillaries) to the 
"root" of the network (e.g., the river outlet or the aorta). 
However, Horton-Strahler ordering cannot be applied to net- 
works with loops because they they create a contradiction in 
the edge ordering in terms of which edge precedes another 
in the hierarchy. Here, we present a generalization of the 
Horton-Strahler order to weighted planar reticular networks, 
where weights are assumed to correlate with the importance 
of network edges, e.g., weights estimated from edge widths 
may correlate to flow capacity. Our method assigns hierar- 
chical levels not only to edges of the network, but also to its 
loops, and classifies the edges into reticular edges, which are 
responsible for loop formation, and tree edges. In addition, 
we perform a detailed and rigorous theoretical analysis of the 
sensitivity of the hierarchical levels to weight perturbations. 
We discuss applications of this generalized Horton-Strahler 
ordering to the study of leaf venation and other biological 
networks. 

Keywords, complex networks, spatial networks, graph the- 
ory, redundancy 



Introduction 

Networks and network theory have been utilized to 
represent and analyze the structure and function of 



*Department of Mathematics, Duke University, Durham, NC 
27708, USA 

t Institute of Science and Technology Austria, 3400 
Klosterneuburg, Austria 

■^Department of Computer Science, Duke University, Durham, 
NC 27708, USA 

§ School of Plant Biology, The University of Western Australia, 
Crawley WA 6009, Australia 

^School of Biology and School of Physics, Georgia Institute of 
Technology, Atlanta, GA 30332, USA 



a myriad of biological systems. These systems span 
scales from cells to ecosystems and include gene regu- 
latory networks [6l |4] , metabolic pathways [27l [22] , dis- 
ease dynamics [38l [36] , food webs [HI [2] , host- parasite 
webs [SH [24], and social interactions [57l [HI [15]. In 
the process, structural archetypes have been identified 
including scale-free behavior, motifs, modularity, the 
emergence of hubs, and small- world structure j57l[7l[52l 
[261 [H im HOI [3]. However, these theories do not typ- 
ically incorporate the spatial constraints that underlie 
the location and connections amongst nodes and edges. 
Indeed, there are many examples of delivery and dis- 
tribution networks where nodes and edges are physical 
structures embedded in space, e.g., leaf venation net- 
works [S31 [TT], cardiovascular networks [351 [55], corti- 
cal networks (SS], root networks [5B1, ant trails [MJ and 
road networks [3S]. Hence, theory is also needed to char- 
acterize biological networks whose structure is strongly 
influenced by physical constraints (for a review, see [8]). 

Although the theory of spatial networks is quite di- 
verse, the theory as applied to resource delivery net- 
works in biology often involves certain simplifying as- 
sumptions. For example, in fractal branching theory, 
a network is seen as a perfectly self-similar structure, 
e.g. a dividing binary tree [46] . A prominent theory of 
metabolic scaling in mammals assumes the cardiovascu- 
lar system is a fractal whose physical dimensions have 
evolved to optimally transport fluid from the aorta to 
capillaries (58] [13] . An extension of this model to the 
above-ground structure of tree branches makes similar 
assumptions |23] . Both models have inspired a wide ar- 
ray of follow-up work with increased recognition that 
the original fractal branching assumption is overly sim- 
plistic [m [44] [50] [5] [31] . For example, in reality, physi- 
cal networks in biology have side branches and are not 
perfectly balanced binary trees [55]. Theories of side- 
branching resource delivery and distribution networks 
have their origins in the study of river networks. In a 
river network, streams merge together to form larger 
streams. However, small streams can merge into larger 
streams of all scales. The topological structure of river 



networks can be analyzed using the so-called Horton- 
Strahler order [28l [51] . This scheme assigns an integer 
number to every branch of the network. The numbers 
represent different levels of the branch hierarchy, with 
larger numbers corresponding to the larger stream seg- 
ments in the network. The Horton-Strahler ordering 
is the basis for the characterization of the statistical 
properties of river networks |20|. including the finding 
that river networks are fractal [17] • Moreover, the side- 
branching statistics first introduced by Tokunaga [SI] 
can be used to characterize universal features of river 
networks and departures thereof [19]. 

Leaf venation networks are a prominent example of a 
physical delivery and distribution network whose struc- 
ture possess numerous side branches. The structure 
of leaf venation networks has broad functional impli- 
cations. For example, leaf vein density is positively 
correlated with photosynthetic rates T^' and also influ- 
ences the extent to which leaves form a hydraulic bot- 
tleneck in whole plants [El [48]. However, many leaves 
of higher plants (notably most leaves of angiosperm 
lineages), have reticulate venation networks, involving 
loops within loops [11 . It has been hypothesized that 
reticulate patterns allow leaves to maintain the supply 
of water and nutrients to and from photosynthetically 
active chloroplasts even when flow through some edges 
in the network is lost [37] |30j |49l [17] due to mechan- 
ical damage or herbivory. Unfortunately, the Horton- 
Strahler ordering scheme developed for the analysis of 
river networks is not directly applicable to reticular net- 
works. The reason is that loops lead to inconsistencies 
in the merging procedure in which a strictly hierarchical 
order is assigned to all edges. 

In this paper we propose a method that generalizes 
the Horton-Strahler order to planar, weighted reticu- 
lar networks. Such networks encompass a large class of 
physical networks, where weights can often be obtained 
by estimating dimensions of edges, such as branch 
widths, or other indicators of cost or importance. While 
coinciding with the Horton-Strahler order for branching 
networks, our method also assigns hierarchical levels to 
the loops of the network. Moreover, it categorizes the 
branches into the ones responsible for the formation of 
loops, and the ones forming the tree structure of the 
network. Edge weights play an important role in our 
algorithm, and we perform a theoretical analysis of pos- 
sible effects of weight perturbations on the hierarchical 
levels. We flnd that the loop hierarchy is more robust to 
measurement error of network edge weights than is the 
tree hierarchy. In the past, comparisons of the statis- 
tical similarity between river networks and leaves have 
been proposed, albeit such comparisons are restricted 
to leaves without loops [43]. Hence, we also discuss ap- 
plications of the current method to the characterization 
and comparison of reticulate leaf venation networks as 



well as obstacles to extending this method to a more 
general class of networks. 

A graph theoretic approach 
to Horton-Strahler ordering of 
rooted trees 

We start by reviewing the algorithm for constructing the 
Horton-Strahler order. For the remainder of the paper, 
we shall adopt the language of graph theory [M] [TO] . 
Note that in graph theory, the "leaves" of the network 
are those vertices which only have a single edge that 
connects to them. In this context, the input to the 
Horton-Strahler ordering algorithm is simply a rooted 
tree, T — {V,E), where V is the set of vertices and E 
is the set of edges. Given such a tree, the algorithm 
assigns a level, A(e), to each edge e e in the following 
way: 

• Assign level 1 to all edges connected to the leaves 
of T. 

• For each vertex having only one incident edge, e, 
with undeflned A(e), let I be the maximal level 
among the other incident edges. If there is a single 
incident edge of level I, then A(e) — I. If there 
are two or more incident edges of level I, then 
A(e) =/ + l. 

The result of this algorithm is illustrated in Fig. [IJa). 
Conventionally in the study of river networks 47 , this 
algorithm can be summarized by a single rule which 
states that the order of a downstream segment is equal 
to 

A = max(Ai,A2) + 6x^^x^ 

where Ai and A2 are the order of the two upstream seg- 
ments that are merging and S is the Kronecker delta. 

It is clear, however, that if the network has loops, as in 
Fig. mb), the algorithm simply cannot proceed because 
there always will be a vertex having more than one in- 
cident edge with an undefined level. Moreover, loops in 
this graph seem to also form a hierarchy. For example, 
the loop outlined in Fig. \V[h) by the red dotted line may 
belong to a higher level than the loop outlined by the 
blue dashed line. It turns out that such a hierarchy can 
be constructed and separated from the tree hierarchy if 
edges have weights and the graph itself is planar. An 
example of such a graph is shown in Fig. [Ijc), where 
the weights represent widths of the branches. 

To guide the reader's intuition, we first provide an 
alternative description of the Horton-Strahler algorithm 
for the case when the tree T is binary and weighted. Let 
w :£'—>■ R be the weight function, that is, w{e) is the 
weight of an edge e £ E. Since the tree is rooted, there 
is a partial order defined on E as follows: ei < 62 if there 
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Fi gure 1 1 Examples of networks with hierarchical structure with 
a common "root" or outlet denoted by the red dot at the bottom 
of each network: (a) Horton-Strahler stream order of branch hier- 
archy in a tree network; (b) Reticular network with possible loop 
hierarchy: the, blue, dashed loop might be less important than 
the red, dotted loop; (c) Reticular network of (b) with weights. 



is a path from the root to ei which passes through 62 
(in other words, 62 is closer to the root than ei). Let 
us assume that the weight function is strictly increasing 
with respect to this order, that is, w{ei) < w{e2) if 
ei < 62- In this case, the Horton-Strahler order can be 
computed using the following procedure: 

• Let C C E he the set of edges incident to leaves 
of T, regarded as a set of disjoint components. For 
each e e C let A(e) — 1. 

• Iterate through (the rest of the) edges in order of 
increasing weight. For each edge e do the following: 

— If e shares a vertex with a single component 
ci G C, then merge ci and e into a new com- 
ponent c, and let A(e) = A(c) = A(ci). 

— If e shares a vertex with two components 
Ci, C2 G C, then merge ci, C2, and e into a new 
component c, and assign levels as follows: 

* If A(ci) A(c2), then A(c) = A(ci) -I- 1. 

* If A(ci) =^ A(c2), then A(c) = 
max{A(ci), A(c2)}. 

* A(e) = A(c). 

Ordering of planar weighted 
graphs 

We now present the algorithm for constructing the gen- 
eralized Horton-Strahler order. Let G = {V, E) be 
a planar graph, not necessarily a tree, and again let 
w : E' — M be a weight function. We shall assume that 
w is injective (i.e., all weights are unique). Otherwise, 
the ties will be resolved arbitrarily. The merging pro- 
cedure for computing the Horton-Strahler order works 
with disjoint components, which, in the language of al- 
gebraic topology, are 0-dimensional homology classes. 




Figure 2: An illustration of the loop merging procedure applied 
to the graph from Fig. [He). Red, dashed edges are the ones 
removed during merging, the corresponding numbers show their 
levels. Levels of faces is encoded by the color: white faces have 
level 1, light blue faces have level 2, and gold faces have level 3. 
Note that /7 is the unbounded face. 



Loops, on the other hand, are 1-dimensional homology 
classes. Hence, we may try to construct a hierarchy by 
merging loops. Notice that the boundary of a face of the 
graph G is a loop, and we can merge two neighboring 
faces by removing a shared edge. Using these two ob- 
servations, we obtain the following merging procedure 
for loops: 

• Sort the edges so that w[ei) < w{e2) < ■ ■ ■ < 
'w(en), where n ~ \E\ is the number of edges. 

• Let A(/) = 1 for each face /. 

• Iterate through ei, . . . , e„ and do the following: 

— If Ci is adjacent to a single face, skip to the 
next edge. 

— If e; is adjacent to two distinct faces /l and 
/r, remove e, from the graph, let /merged = 
/l U fji, and assign the levels as follows: 

* If A(/l) = Xifn) then \{f,nerged) = 

Hh) + 1. 

* If A(/l) ^ XUr) then X{f,nerged) = 

max{A(/i),A(/«)}. 

* A(e,) = min{A(/i),A(/fl)}. 

A step-by-step illustration of loop merging applied to 
the tree in Fig. [ijc) is shown in Fig. [2] Notice that 
this procedure builds a rooted binary tree, where leaves 
correspond to the faces of G, and the rest of the vertices 
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correspond to unions of these faces. The assignment of 
levels in this tree follows the original Horton-Strahler 
algorithm. It is also useful to remember that faces of G 
are vertices of its dual graph, G* , and merging faces of 
G can be thought of as adding an edge to G*. Hence, 
the two merging procedures that we described are, in 
some sense, dual. We shall refer to the binary tree of 
faces as the co-tree of G, and denote it by T*{G). 

The construction of T*{G) removes edges from G 
which are responsible for the existence of loops. We 
shall call such edges reticular. Assignment of levels for 
such edges is based on the assumption that a merger 
should not be more significant than any of the merging 
elements. Notice that after removing reticular edges 
from G we have a spanning tree of G, which we denote 
by T{G). This tree captures the tree- like structure of 
the original network, and we can assign hierarchical lev- 
els to its edges using the original Horton-Strahler algo- 
rithm. We only need to determine which vertex should 
be the root, and we do this by finding the vertex with 
a single incident weight of maximum weight. Hence, we 
augment the procedure for constructing the loop hier- 
archy by the following statement: 

• Apply the Horton-Strahler ordering to the remain- 
der of the graph (which is a rooted tree). 

The result of the complete algorithm applied to the tree 
in Fig. [TJc) is provided in Fig. [3l 

The algorithm produces three types of output. First, 
it provides a unique set of orders to those edges involved 
in the non-reticulate component of the network (Figure 
3 - left panel). Second, it provides a unique set of or- 
ders to those edges involved in the formation of loops 
(Figure 3 - right panel). Further, one can also calculate 
the side-branching statistics associated with both or- 
derings. The side-branching statistics, i.e., "Tokunaga" 
statistics [541 , for a conventional non-loopy tree are sum- 
marized by the numbers Nij which are the number of 
edges of level j that join with edges of level i. Because 
of the ordering process, these statistics are evaluated for 
i > j. These numbers can also be divided by the num- 
ber of absorbing edges, i.e., the total number of edges 
of level i to yield an average number of side-branches 
per segment. Here, the algorithm produces two sets 
of Tokunaga statistics, the numbers Nij for the side- 
branching of tree edges (Figure 3 - left panel) and Mij 
for the side-branching of reticulate edges (Figure 3 - 
right panel). 

Sensitivity of planar network or- 
dering to weight perturbations 

Clearly, edge weights play an important role in the con- 
struction of both loop and tree hierarchies. Unfortu- 




A'2i = 2 A'j2 = jVj, = 1 M21 = 1 M,2 =1 A/ji = 

Figure 3: Hierarchical levels assigned to the loops and branches 
of the network from Fig. [TJc). Edge levels are shown on the left, 
where black edges have order 1, light blue edges have order 2 and 
gold edges have order 3; reticular edges are dashed. Face levels are 
shown in the co-tree on the right, where white nodes have order 1, 
light blue nodes have order 2 and gold nodes have order 3. Leaves 
of the co-tree are labeled by the corresponding faces while other 
nodes are labeled by the reticular edges causing the merger of the 
two child nodes. Numbers A''^^ are Tokunaga statistics for the 
spanning tree and indicate the number of edges of level j joining 
with edges of level i |54j . Similarly, Mij are Tokunaga statistics 
for the reticulate co-tree and indicate the number of edges and 
faces of level j merging with edges of level i. For both M and A'^, 
statistics are only collected when i > j. 



nately, weight estimation done in practice is often im- 
precise, so the order in which the algorithm iterates 
through the edges may be perturbed. In this section 
we investigate how such a perturbation afi'ects the loop 
and tree hierarchies. 

We start by considering the worst possible change in 
the hierarchical levels of loops. Notice that the highest 
level in the hierarchy of loops can be as low as 2. This 
happens when the first reticular edge creates a level 2 
face and every other reticular edge merges a level 1 face 
with the only level 2 face (see Fig. |H[a) ). On the other 
hand, the highest level in the loop hierarchy can be as 
high as log(m), where m is the number of faces. This 
happens when level 1 faces are merged only with level 
1 faces until only faces of level 2 are left, then level 2 
faces are merged with level 2 faces until only faces of 
level 3 are left, and so on (see Fig. SKb) ). It is clear 
from the example in Fig. |4]that there is a permutation of 
edges that can change the loop hierarchy from one of the 
extreme cases to the other. However, in practice such a 
permutation would generally result in from a significant 
perturbation in weights. For small perturbations, it is 
more likely that only a few transpositions of edges will 
occur. 

Let ei, . . . , e„ be the order of edges with respect to 
their weights. We shall now analyze how the structure 
of T{G) and T*{G) changes when a single transposi- 
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Figure 4: Example of the two extreme cases of the loop hier- 
archy. The network has m faces, where m = 2^ for some integer 
fe > 0. (m — 1) of these faces are adjacent squares and the other 
one is the unbounded face. Vertical edges are removed before hor- 
izontal edges as follows; (a) The edges are removed sequentially 
from left to right. The corresponding co-tree has the shape of a 
"comb" and the maximal hierarchical level is 2; (b) The edges are 
removed from left to right skipping every second edge. The pro- 
cess is repeated until all vertical edges except the rightmost one 
are removed. The corresponding co-tree has the height logira) = k 
which is the maximal hierarchical level. 



tion occurs, that is, when the order of and e^+i is 
swapped. First, we notice that there will be no changes 
to the structure of the co-tree or the spanning tree if 
Ci and Ci+i are both tree edges, or if is a tree edge 
and Ci+i is a reticular edge. Hence, there are two cases 
to consider: when both and e^+i are reticular, and 
when Ci is a reticular edge and e^+i is a tree edge. In 
the former case, we can regard reticular edges as edges 
of the co-tree. We see then that swapping the two edges 
may shift a subtree of the co-tree only one level up or 
down. Therefore, it is reasonable to expect that hierar- 
chical levels of loops will change at most by one. The 
case of a reticular edge and a tree edge is more com- 
plicated. Such a transposition may lead to detaching a 
subtree of the remaining spanning tree and attaching it 
at a different place. This may have a drastic effect on 
the tree hierarchy. A detailed analysis of the two cases 
justifying the above conclusions is present below. 



Case 1. Ci and e^+i are both reticular. Only the co- 
tree can be affected in this case. Let //^, and fHi, 
ft+i be the faces merged by removing ei and ei+i, re- 



spectively. Also, let /, = U/« and = ff+.UfH^. 
Notice that if fi ^ //^^i and ^ //^j, then fi is not 
a child of fi^i in T*{G)), and there will be no changes 
to the structure of the co-tree. Suppose that fi = f^_^-^ 
(the case when fi = f^i follows the same argument). 
Then e^+i is adjacent to either //" or //*-; let us assume 
it's //^. Removing e^+i before leads to merging f^^ 
with //^ first, and then merging the resulting face with 
f^. The corresponding change in the tree structure, 
shown in Fig O is a single rotation around /i+i- Pos- 
sible changes in the levels of the nodes involved in the 
rotation are also shown in Fig. [5l We can see that these 
levels can change at most by one. However, in the worst 
case the change in levels may propagate up T*{G) all 
the way to the root. 



Case 2. is a reticular edge and e^+i is a tree edge. 
Let //^ and //^ be the two faces merged by removing e^. 
Notice that there will be no changes in the structure of 
T(G) or T*{G) if e^+i is not adjacent to both //- and 
f^. So, let Ci+i be adjacent to f^ and f^. Then re- 
moving e^+i before merges the same f^ and //^, so 
no changes to the structure of the co-tree happen. How- 
ever, e^+i turns into a reticular edges, and e,; becomes a 
tree edge. Consequently, the structure of the spanning 
tree changes. Let E^'^ be the set of edges incident to 
both fl" and //^, and let T^'^ be the tree formed by 
the edges in E'"^^ and the edges connected to iJ^-^ and 
having only or as an adjacent face (see Fig. [6]). 
Removing e^+i and splits T^'^ into three trees, T^, 
T^, and T^^, such that and are connected to 
the boundary of //^ U /f , and T^^ is not (Fig. [6]). If 
is removed before e^+i, then T*^ is connected to by 
Bi+i. However, if the transposition happens and e^+i 
is removed before ei, then T^^ is connected to by 
ei (Fig. [B]). To understand the effect of such a change 
on hierarchical levels, we first assume that T^' does not 
contain the root of T{G). Let be the vertex incident 
to ei+i and T^, and let be the vertex incident to 
ei and T^. Also, let A^+i = A(ei+i), where Ci+i is re- 
garded as an edge in T^^ U Ci+i rooted at v^, and let 
\i — A(ei), where is regarded as an edge in T*^ U 
rooted at . Denote by the edge of T{G) which is 
next to ei+i in the path from ei+i to the root of T{G), 
and by the edge of T{G) which is next to ei in the 
path from to the root of T{G). Then we can see that 
removing e^+i before can decrease the level of by 
at most max{l, A^+i — 1}. At the same time, the level of 
can increase by at most max{l, A^ — 1}. In the worst 
case, these changes can propagate up T{G) all the way 
to the root. The case when the root of T(G) belongs 
to T*^ can lead to more drastic changes. In this case, 
removing Ci+i before leads to recomputing levels of 
all edges in T{G) — T^'^ by changing the root from 
to . Again, this change can then propagate further 
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Figure 5: Example of the effect of a single transposition of two reticular edges: (a) the part of the network containing the two edges 
being transposed and the effect of the transposition on the structure of the co-tree; (b) possible level changes caused by the transposition. 




Figure 6: Example of the effect of a single transposition of a reticular edge and a tree edge: (a) The part of the network containing the 
two edges being transposed. The brown, blue, and green triangles (and edges) denote the subtrees adjacent to the edges, (b) The effect 
of the transposition on the structure of the spanning tree and its hierarchical levels. 



to the root of T(G). 

Discussion 

We have shown that the hierarchy of loops often ob- 
served in reticular physical networks can be defined ex- 
plicitly using a generalization of the Horton-Strahler or- 
der. To obtain such a generalization we regard the net- 
work as a weighted graph, with weights corresponding 
to the widths of the network branches. Noticing that 
the Horton-Strahler order can be computed by analyz- 
ing how specific disjoint components (sub-networks) of 
a (non-reticular) network are merged as the edges are 
added in the order of increasing weight, we show that the 
hierarchical order of loops in a weighted planar graph 
can then be computed by analyzing how the faces of the 
graph are merged as we remove the edges in the order 
of increasing weight. This approach naturally classifies 
graph edges into reticular edges, which are responsible 
for loop formation, and tree edges, which constitute a 
spanning tree of the graph. Hence, both the loop and 
the tree hierarchies can be computed. 

Being able to compute hierarchical levels for loops 
creates new possibilities for analyzing the structure of 
reticular networks. By means of analogy, river networks 



can be compared by representing their connectivity in 
terms of side-branching statistics |54] . These statistics 
depict the ways in which smaller streams connect to 
larger streams at all scales of the network [19] . A sim- 
ilar procedure could be applied to leaf networks. For 
example, the current algorithm decomposes reticulate 
networks into a binary-tree for the loop hierarchy and 
a separate binary tree for the tree hierarchy. Both net- 
works have associated Horton-Strahler orders and there- 
fore their structure can be estimated using Tokunaga 
statistics. Recent innovations in software now permit 
semi-automated extraction of the dimension and con- 
nectivity of entire leaf vein networks and the areoles 
that veins surround |45| . Hence, greater quantification 
of leaf vein networks from across a wide range of biolog- 
ical diversity will soon be available for which to analyze 
leaf development, variation across environmental gradi- 
ents and in paleobotanical studies. Current attempts 
to compare reticulate structure have largely focused on 
the density of areoles (i.e., network faces) as a proxy for 
the "loopiness" of the network [9]. The current study 
will provide additional metrics to compare the detailed 
branching structure of reticulate networks. 

An important caveat to keep in mind when comparing 
reticulate network structure is that estimating weights 
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in physical networks is by no means a trivial problem. 
Therefore, we have performed a theoretical analysis of 
possible changes in the loop and tree hierarchies due to 
perturbations in edge weights. We have shown that the 
worst possible change in the loop hierarchy is attain- 
able, but requires a significant perturbation of weights. 
Taking into account that small perturbations are likely 
to cause only a few transpositions in the order in which 
the edges are removed, we have shown that a single such 
transposition can change the hierarchical levels of loops 
at most by one. We have also shown that the change 
in the hierarchical levels of the remaining spanning tree 
can be arbitrarily large even when a single transposi- 
tion is performed. It is important to note that in either 
case the change does not happen for every transposition. 
Rather, the transposed edges have to satisfy a partic- 
ular condition, which may happen rarely in practice. 
The latter claim is supported by the numerous success- 
ful applications of the Horton-Strahler order. While the 
method itself does not depend on any weights, the con- 
nectivity of the network is obtained by analyzing digital 
elevation map which contain noise [531 IH] • In particu- 
lar, the difference between the correct and the computed 
connectivity may be exactly the same as the difference 
in the connectivity of our spanning tree caused by trans- 
posing two edges. Hence, the resulting hierarchy may be 
drastically different from the correct one. Nevertheless, 
the Horton-Strahler order has been successfully used for 
over five decades despite the potential instability iden- 
tified here UHllSIl [53113 US- We suggest that empirical 
characterizations of reticulate planar networks include 
randomization analysis on edge weights to identify the 
robustness of claims regarding statistical structure of 
side-branching of the tree and co-tree. 



Many biological and physical systems are represented 
by non-planar physical networks ,25, 8 and comput- 
ing hierarchical levels of loops in such networks is still 
an open question. While our method can be applied 
to obtain the tree hierarchy of such networks, the loop 
hierarchy cannot be computed in this case because the 
algorithm relies on the fact that any loop in a planar net- 
work corresponds to a union of faces. In the mathemat- 
ical language, (boundaries of) faces of a planar graph 
form a canonical basis for loops (1-dimensional homol- 
ogy classes). Such a canonical basis is not present in 
non-planar graphs. It is not clear at this point how 
to handle the non-planar case. Perhaps a method for 
computing loop hierarchies which is not based on the 
widths of the network branches could provide an an- 
swer. We hope that our approach of using algebraic 
topology language to deal with nodes and loops of net- 
works will prove useful in developing such a method and 
complement other approaches. 
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